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Elastic scattering of photons in a Lorentz-scalar potential via virtual spin-zero 
particle-antiparticle pairs ("Delbriick scattering") is considered. An analytic expres- 
sion for the Delbriick amplitude is found exactly in case of an oscillator potential. 



_C ■ General properties of the amplitude and its asymptotics are discussed. 
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1. Elastic scattering of photons in the Coulomb field of nuclei via virtual electron- 
positron pairs (Delbriick scattering |I[) attracted considerable interest for a long time. 
That is motivated by two reasons: (i) The Delbriick scattering is one of a few nonlinear 
quantum-electrodynamic processes which can be precisely tested by experiment || 0, f|, ^]. 
(ii) In order to extract an information on nuclear structure from differential cross sections 
of photon scattering on nuclei, a precise knowledge of the Delbriick amplitude may be 
required because of its interference with the nuclear amplitude. In some cases, the Delbriick 
scattering considerably modifies the differential cross section of 7A-scattering || 0, ||, || . 

Theoretical investigations of the Delbriick scattering have a long history and many 
papers are devoted to this subject. Now the Delbriick amplitude is studied in detail in 
some approximations: (1) In the lowest-order Born approximation with respect to the 
parameter Za (here Z\e\ is the charge of the nucleus, a = e 2 ~ 1/137 is the fine-structure 
constant), results were obtained for an arbitrary momentum transfer q; these results are 
surveyed in detail in ||. (2) For the case of high energies (uo ^> m e , m e being the electron 



mass) and small scattering angles (q <C lu), the Delbriick amplitude was obtained in |T0, 
|TT1 , |12"1 , |13| , [HH to all orders in Za. It was found that the Coulomb corrections at Za ~ 1 
drastically change the amplitude as compared to the Born approximation. (3) At high 
energies and momentum transfers (u, q ^> m e ), the amplitude was also found to all orders 
in Za JT5], [T5], |T7[]. In this case the Delbriick amplitude has a scaling behavior [I8| and 



becomes inversely proportional to the photon energy. It has been shown that the Coulomb 
corrections essentially decrease the amplitude at high momentum transfer as well. Recently 
some general exact expressions for the Delbriick amplitude were derived in |19|], although 



numerical results were not presented. Some numerical results for the Delbriick amplitude 
obtained in the lowest-order Born approximation, in the "high-energy and small-angle" 
approximation, and in the "high energy and momentum transfer" approximation can be 



found in [ED]. 



In the present paper we consider the Delbriick scattering in the scalar QED in the field 
of an oscillator Lorentz-scalar potential, or, in other words, for a relativistic oscillator. The 
relativistic oscillator model, being a nice theoretical laboratory, has also a few realistic 
applications. For instance, it can be successfully used |2l[ to describe interactions of col- 
lective modes of a nucleus in the region of giant resonances with the electromagnetic field. 
It takes the advantage of being automatically consistent with microcausality, analyticity, 
and dispersion relations. The Delbriick amplitude is a quantum correction to a classical 
part of the photon scattering by the oscillator and it is interesting to evaluate this correc- 
tion explicitly to understand its significance for the physics of photon-nucleus scattering. 
Besides, the Delbriick amplitude for the relativistic oscillator provides an example of exact 
calculation of this quantity in external potentials. 

2. To be specific, we consider spin-zero particles described by the Klein-Gordon equa- 
tion with an oscillator Lorentz-scalar potential. We start with the Lagrangian for the 
charged quantum field <fi in the external electromagnetic potential A^ : 

L[4>,A](x) = \d^(x)+ieA ll (x)4>(x)\ 2 -U(r)\^(x)\\ U(r) = fi 2 + 7 V. (1) 

In the Furry representation, the Feynman rules corresponding to this Lagrangian are found 
from the usual Feynman rules for spin-zero particles by replacing plane waves with the 



normalized solutions (r\n) = <f) n (r) of the Klein-Gordon equation in the potential U(r) : 

p 2 Ur)+l'r 2 Mr) = (E 2 n - l x 2 )Mr), P = -«V (2) 

(for brevity, we do not indicate here quantum numbers related to the angular momentum). 
The equation (Q) coincides with that for a nonrelativistic oscillator and has the spectrum 
E = ±E„, where 



J n- 



E n = y/fj? + (2n + 3) 7 2 , n > 0. (3) 

The parameter 7 4 determines a slope of the potential U(r) in ([I]) and must be positive 
to result in a stable vacuum and bound states. It also determines a range of the ground 
state, 0o( r ) — 7 3 ' 2 7r _3 ' 4 exp(— |7 2 r 2 ). The parameter /i 2 may be negative. However, we 
require /i 2 > — 37 2 to have a positive energy gap 2E Q between levels of positive and negative 
energies and hence a stable vacuum. In the nonrelativistic limit, 7 -C /i, the parameter \x 
becomes the mass m of the particle and the oscillator parameter 7 determines the oscillator 
frequency u Q = E n+ i — E n =const, 

H — > m, 7 2 — * mujQ. (4) 

The Delbriick amplitude for the spin- zero particle is described by two Feynman diagrams 
shown in Fig. 1. In this figure, double lines represent the Green function for the Klein- 
Gordon equation in the external potential U(r). The first diagram, Fig. la, describes the 
so-called seagull contribution to the photon scattering amplitude which is equal to 



/OO g^g p I 

— / d 3 rG(r,r\e) exp(iqr) = — a(ee') V] — -(n\ exp(iqrr)|n). (5) 
-00 2n J „ E„ 



Here e (k) and e' (fe ) are polarizations (momenta) of the incoming and outgoing photons, 
respectively, and q = k — k' ; we use the radiative gauge, ek = e'k' = 0. G(r, r'\e) is the 
Green function in the potential U(r): 

. . (r\r)\ (ti 7* ; r 00 r n 

G ( r > r '^ = E 2 La \ n = - % \ ds ^ ex P H" 2 -P 2 - U ^ + ? °) l r ')' ( 6 ) 

n e ~ E n + lQ Jo L J 

s being a proper time. The amplitude we use is normalized as to give, being squared, 
exactly the differential cross section of photon scattering. 
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Figure 1: Diagrams of the Delbriick scattering: (a) Seagull amplitude So- (b) Resonance 
amplitude Rd- 



Using a known expression for the time-dependent Green function of a nonrelativistic 
particle in an oscillator potential (see e.g. |22[), we obtain an explicit form for the Green 
function (|6f): 



G(r,r'\e) 



-Z7 



is 



ds 



exp 



(27risin2s) 3 / 2 
e 2 - /i 2 + zO) + 



x 



27 



r' ) cos 2s 



2rr' 



v 7^ ' 2 sin 2s 

Substituting (^) to (pj) and taking elementary integrals with respect to e and r, we find: 



(7) 
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(8) 



In case of the oscillator potential U(r), the Klein-Gordon operator has only a discrete 
spectrum (|3|) and therefore the Green function has only simple poles at real e = ±(E n — iO) 
displaced by iO in accordance with the Feynman rules. Using this analytic property, we 
can deform the integration path over e in (|J), e — > ie, to make it finally coincident with the 
imaginary axis. As a result, we get a possibility to rotate the contour of integration over s 
in ©-(H), s — > —is, to make all the integrals well-convergent. 

Note that an ultraviolet divergence at s — > is not present in the amplitude Sd provided 
q 7^ 0, so that a regularisation is not necessary. The same is valid for the second diagram, 
Fig. lb, which gives the so-called resonance contribution: 
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ie''V'G(r',r\e + uS) exp(ikr — ik'r') 
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(9) 
(10) 



It has a series of poles at uj = ±(E n + E n >) which are related with the pair photoproduction 
from vacuum when the pair components are captured to the levels n, n! of the discrete 
spectrum. The pole at u = 2E is nevertheless absent because matrix elements in (|10|) 
vanish for the s-wave states n = n' = 0. Using (|7|), we easily calculate the integrals (§) 
with respect to the variables e, r and r' (here the substitutions r = r\ + T2 and r' = ri~r 2 
are helpful). It is also convenient to replace Si = s(l + x)/2 and s 2 = s(l — x)/2, where Si 
and S2 are proper times in the integral representation of the Green functions in (|9|). Then 
we get 
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where 



F = -s 



fi 2 -i0 
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sin [sx) 
2j 2 sin s 



[ek')(e'k) — icos(sx)(ee') 
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q z -2uo z . . uj z 
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(11) 



(12) 



Integrating the term in (|TT|) proportional to (ee') by parts with respect to the variable x, 
we find that the term outside the integral cancels the contribution Sd of the first diagram, 



so that the Delbriick scattering amplitude Tjj = So + Re> finally reads: 

ae^^uj 2 f°° r 1 <Jssin(sx) , ... ... 4F 

Td = —-—~z~j=I ds dx e isms (ee ) — sm(sx)(ss ) e , (13) 

167 ci v 7r "'o "'o sin s 

where (ss 1 ) = (fe x e)(fc x e') describes a magnetic response. We see that a contribution 
like >Sd which depends on q and is independent of uj vanishes in the total amplitude Tjj, as 
it must be according to general consequences of gauge invariance and related low-energy 
theorems. Eq. ( fl3T) is our main result. 



We may obtain another form of Eq. (13) by deforming the integration path over e in 

(0) to transform the integral to j_? _. de. Such a deformation is always possible at low 

energies uj because the chains of singularities of two Green functions in @ lying below and 
above the real axis, i.e. at e = E n — uj — iO and at e = — E n + iO, respectively, do not pinch 
the integration path provided uj < 2E . Then we can rotate si — > — is\, s 2 — > — is 2 and 
get well-convergent integrals. Respectively, we may rotate the integration path in flT3]), 
s —>■ —is, and arrive at a real integral which is well suitable for calculating Tp at low 
energies. However, such an integral turns out to be divergent at uj > E Q + Ei, i.e. above 
the nearest pole of the amplitude T D , see Eq. fllPP- 



3. Eq. ( |I3|) can be further simplified in case of a low oscillator frequency and momentum 
transfer, 7 2 <C fJ 2 and q <^C /x. Then the contribution of the region s > 1 to the integral in 



(|13|) is exponentially suppressed, as is seen after the rotation of the contour, whereas the 



contribution from the region sC 1 reads 

ae in / 4 uj 2 r , , r°° ds r 1 
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i6 7 3^ [(ee/) " (ss/)] r^/o x2 ^ x 

,2 „-n „2 /i ^2N2 



cxp 



. fi^-iO . q z . {1-x- j 2 3 
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7 2 4 7 2 s 487 2 



(14) 



and is saturated by 



I g/i x , if fJ,q>T- 



Respectively, at "low" energies uj 2 <^C 7 2 s^| which may nevertheless be high in comparison 
with /i, the integral ( ]T4| ) reduces to the modified Bessel function of the third kind K 3 / 2 and 
hence is an elementary function: 

(16) 

When q — > 0, the Delbriick amplitude becomes proportional to an effective volume of the 
vacuum, V ~ g -3 , probed by the photons. At high energies uj 2 3> 7 2 s^, the integral (O) 
is saturated by x close to 1 and the answer is given by the modified Bessel function K 3 : 




D 



■a.\x 
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Tq 




±iujV3^K 3 (q\[(ee')-(88% (17) 



where the sign of plus or minus appears for Imuj > or Imw < 0, respectively. Since 
K 3 (z) — > 8z~ 3 when z — > 0, the last expression seems to have a too strong singularity 



~ q 6 ~ V 2 when q — > 0. However, when q is very small, the assumption uj 2 ^> 7 2 s ef | 
becomes violated and the Delbriick amplitude approaches to the regime of Eq. (p~6|) . 



The singular behavior of the Delbriick amplitude at small q, Eq. fllq) , can be understood 
in the following way. If the external potential U is equal to zero, the photon scattering 
amplitude due to the diagrams shown in Fig.l reads [p3| : 



• A -^ e "* eV (K k v ~ &9ia>) log — 2 : 



T = (27r) 3 5 3 (fc - fe'^e'V^ - k 2 g, u ) log — 2 , (18) 



where A is an ultraviolet cut off and m is the mass of bare scalar particles, and is pro- 
portional to w 2 [(ee') — (ss')]. That means that the vacuum has an electric and magnetic 
polarizability per unit volume, 

A 2 

vac.O vac.O " i **■ /in\ 

Xe =~Xm =^log^|- (19) 

These universal vacuum susceptibilities are absorbed by renormalisation of electromag- 
netic fields and charges and are not observable. However, in the presence of an almost 
uniform scalar potential U(r) which shifts the mass of the particles m 2 , — > m 2 ff = U(r), the 
polarizabilities get a finite meaningful piece, 



a U{r) 

48^ iog "^g" 



xT(r) = -XlT(r) = ~— 2 log -V. (20) 



The last formula is valid when the potential U(r) is constant at distances ~ m e g which are 
characteristic for creating those vacuum polarizabilities: 



mj < Ar = U(r) 



dU{r) 
dr 



(21) 



For the potential U(r) = /i 2 + 7 4 r 2 , it means 

J r > 7 _1 , if 7 2 r >// 



r < ^i 3 7 4 , if 7 2 r < // 



In the case 7 <C //, both the above regions overlap and the formula (|20f) is valid everywhere. 
Then the scattering amplitude of low-energy photons is equal to 

T D = uj 2 a E (q)(ee') + w 2 a M (9)(^') (23) 

with 

d f U(t) 

a E (q) = -a M {q) = -7—^ / exp(iqr) log — r d 3 r . (24) 

487T 7 tuq 

Taking this elementary integral, we find 

a E (q) = -a M (q) = C5 3 (q) + -^- 3 (l + a)e- a } a=^, (25) 

Ylq 6 7 Z 

in complete accordance with (0). Here C is an infinite constant due to the polarizability 
of the vacuum in the whole space; it is infinite because the potential U(r) — > 00 when 

r — > 00. 

6 



The approximation of a uniform potential (|20f) turns out to be inapplicable at high 
energies when the virtual particles of the mass m c s produced by the photon propagate to 
a distance ~ he/ AE ~ u/m 2 s which is large in comparison with the scale Ar of variation 
of the potential U(r). This just happens when uo 2 ^> 7 2 s^f and the amplitude To becomes 
predominantly imaginary, see Eq. fllT" 



Keeping in (f[3|) terms of the next order in s 2 , we can find a correction to the expansion 
of Td in powers of the momentum transfer. At "low" energies it reads 



auj 2 e c 



T D ^Eq. (1Q) + ———^(19 + 2a)(ss') - (17 + 2a)(ee')] . (26) 

It determines an asymptotics of the helicity-non-flip amplitude T^ + because the piece 
(O) contributes to only the helicity-flip amplitude Tp, as easily seen from the relation 
(ee') = ±(ss') = (1 ± cos#)/2 for the helicity-non-flip and helicity-flip case, respectively, 
9 being the scattering angle. 

The singular behavior of Td at q — > disappears when the potential U(r) has a finite 
range. For example, in case of a cut oscillator potential, 



^--Aexpf-^" 



log^ = -Aexp(-^H, (27) 



U(0) \m 

polarizabilities (|24D are finite: 



where A = log — -— = O ( — ) ^C 1 and ujq is the frequency of small oscillations, the 



aA 5 / 2 (_Af 



«b(s) = -®M{q) = AO 7=. , 3 exp ( -^-^ . (21 



At energies uo = 0(ujq) the corresponding Delbriick amplitude (|23|) is less by the fac- 
tor of O((u>o/mo) 3 / 2 ) than the ordinary nonrelativistic photon scattering amplitude Tq ~ 



auj 2 



m (u)l - u 2 ) 



(ee 1 ) by a particle bound at the ground state. 



4. In conclusion, in the present paper we calculated for the first time the amplitude 
of the Delbriick scattering in the scalar QED in case of a scalar oscillator potential and, 
at low oscillator frequency, investigated its asymptotic behavior at low and high energies 
and low momentum transfer. A close relation of the Delbriick scattering with the vacuum 
polarisation was demonstrated. 
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